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Abstract 

We analyze the difFerential equations produced by the method of creative telescoping applied 
to a hyperexponential term in two variables. We show that equations of low order have high 
degree, and that higher order equations have lower degree. More precisely, we derive degree 
bounding formulas which allow to estimate the degree of the output equations from creative 
telescoping as a function of the order. As an application, we show how the knowledge of these 
formulas can be used to improve, at least in principle, the performance of creative telescoping 
implementations, and we deduce bounds on the asymptotic complexity of creative telescoping 
for hyperexponential terms. 
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1. Introduction 

Creative telescoping is a technique for computing difFerential or difference equations 
satisfied by a given definite s um or integral. The technique became widely known through 



the work of lZeilbergeii (|1991[ ). who first observed that creative telescoping in combination 
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Fig. 1. Sizes {r, d) of creative telescoping relations for the integral of a certain rational function 



with Gosper's algorithm (jGosperl . 119781 ) for indefinite hypergeometric summation leads 
to a complete algorithm for computing recurrence equations of definite hypergeom et- 
ric sums. This algorithm is now known as Zeilberger's algorithm (|Zeilbergeii [l990[) . In 
its original version, it accepts as input a bivariate proper hypergeometric term f[n,k) 
and returns as output a linear recurrence equation with polynomial coefficients satisfied 
by the s um F{n) = J^t-^ fi^jk). An an alogous algorithm for definite integration was 



given by lAlmkvist and Zeilberger ( 1990l ) . This algorithm accepts as input a bivariate 
hyperexponential term /(x, y) and returns as output a linear differential equation with 
polynomial coefficients satisfied by the integral F{x) — / (a;, y)dy. A summary of the 
method of creative telescoping for this case is given in Section 2 below. For fu r ther de- 
tails, variations, and gen e ralizations, consult for instance IPe tkovsek ct al. ( 1997 ) , ChvzakI 
(|2000t ) , ISchiieideil |2005l).IChvzaket all 1200^ iKauers and Paulc (20111). For im p lemen - 
tations, see |Paule and SchornI (ll99^7lOffizay (|l998[) . iKoepj (|l998l) . ISchneideil (|2004[ ). 
Abramov et al] (|2004[ ) . iKoutschanI (|2009l l2010l ). etc. 

The equations which can be found via creative telescoping have a certain order r and 
polynomial coefficients of a certain degree d. But for a fixed integration problem, r and d 
arc not uniquely determined. Instead, there arc infinitely many points (r, d) e such 
that creative telescoping can find an equation of order r and degree d. These points form a 
region which is specific to the integration problem at hand. Figure 1 shows an example for 
such a region. Every point (r, d) in the gray region corresponds to a differential equation 
of order r and degree d which creative telescoping can find for integrating the rational 
function 



(Sx'^y^ + 9x^y + 9x'^ + lOxy^ + 3xy + Ax + ij j/ (ix^y^ + 9x^y^ + x^y + 3x 
+ Tx^y^ + Sx'^y'^ 



5x 



IQxy^ + 3xy 

8xy^ + lOxy^ 



lOxy 



by" + lOy^ 



by 



The picture indicates that low order equations have high degree, and that the degree 
decreases with increasing order. But what exactly is the shape of the gray region? And 
where does it come from? And how can it be exploited? These are the questions we 
address in this article. 
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How can it be exploited? There are two main reasons why the shape of the gray 
region is of interest. First, because it can be used to estimate the size of the output 
equations, and hence to derive bounds on the computational cost of computing them. 
Secondly, because it can be used to design more efficient algorithms by recognizing that 
some of the equations arc cheaper than others. 

An analysis of this kind was first undertaken by iBostan et al. I (|2007l ). They studied 



the problem of computing differential equations satisfied by a given algebraic function 
and found a similar phenomenon: low order equations have high degree and vice versa. 
Among other things, they found that an algebraic function with a minimal polynomial of 
degree n satisfies a differential equation of order at most n with polynomial coefficients of 
degree 0(n'^), but also a differential equation of order 6n whose coefficients have degree 
only O(n^). Their message is that trading order for degree can pay off. 



The same phenomenon applies to creative telescoping, as was shown by IBostan et al 



(I2OIOI) For the case of integrating rational functions. The results in the present article 



extend this work in two directions: First in that we consider the larger input class of 
hyperexponential terms, and second in that we give not only isolated degree estimates 
for some specific choices of r, but a curve which passes along the boundary of the gray 
region and thus establishes a degree estimate as a function of the order r. 

Where does it come from? The standard argument for proving the existence 
of creative telescoping relations rests on the fact that linear systems of equations with 
more variables than equations must have a nontrivial solution. Every creative telescoping 
relation can be viewed as a solution of a certain linear system of equations which can be 
constructed from the data given in the input. There is some freedom in how to construct 
these systems, and it turns out that this freedom can be used for making the number of 
variables exceed the number of equations, and thus to enforce the existence of a nontrivial 
solution. 

This reasoning not only implies the existence of equations and the termination of the 
algorithm which searches for them, but it also implies bounds on the output size and 
on the computational cost of the algorithm. But in order to obtain good bounds, the 
freedom in setting up the linear systems must be used carefully. For a good bound, we 
not only want that the number of variables exceeds the number of equations, but we also 
want this to happen already for a reasonably small system. The shape of the gray region 
or iginates from the smalle st systems which have solutions. 

Verbaete"nl ( 1974 . 1976 ) introduced a technique which helps in keeping the size of the 



systems small. The idea is to saturate the linear systems by introducing additional vari- 
ables in a way that avoids increasing the number of equations. We will make use of this 
idea in Section 3 where we propose a design for a parameterized family of linear sys- 
tems whose solutions give rise to creative telescoping relations. Unfortunately, it requires 
some quite lengthy and technical calculations to translate this particular design into an 
inequality condition which rephrases the condition "number of variables > number of 
equations" in precise terms. However, as a reward wc obtain a good approximation to 
the gray region as the solution of this inequality. 

What is the exact shape? We don't know. All we can offer are some rational 
functions which describe the boundary of the region of all (r, d) where the ansatz de- 
scribed in Section 3 has a solution (Theorem 14) . The graphs of these rational functions 
are curves which pass approximately along the boundary of the gray region. 
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By construction, for all integer points (r, d) above these graphs we can guarantee 
the existence of a creative telescoping relation of order r with polynomial coefficients 
of degree d. But we have no proof that our curves are best possible. Experiments have 
shown that at least in some cases, our curve describes the boundary of the gray region 
exactly, or within a negligible error. In other cases, there remains a significant portion of 
the gray region below our curve when r is large. 

In cases where the curve from Theorem 14 is tight, we can compute the points (r, d) 
for which certain interesting measures (such as computing time, output size, . . . ) are 
minimized, as shown in Section 5. Even when the curve is not tight, these calculations 
still give rise to new asymptotic bounds (including the multiplicative constants) of the 
corresponding complexities. We expect that this data will be valuable for constructing 
the next generation of symbolic integration software. 



2. Creative Telescoping for Hyperexponential Terms 

We consider in this article only hyperexponential terms as integrands. Throughout 
the article, K is a field of characteristic 0, and ]K{x,y) is the field of bivariate rational 
functions in x and y over K. Let and Dy denote the derivations on ]K.{x,y) such 
that DxC ~ DyC = for all c e K, and D^x = 1, D^y = 0, DyX = 0, Dyy = 1. 
One can sec that and Dy commute with each other on ]K{x, y). Wc say that a field IE 
containing K(a;, y) is a differential field extension of K(a;, y) if the derivations and Dy 
are extended to derivations on IE and those extended derivations, still denoted by D^ 
and Dy, commute with each other on E. 

Definition 1. An element h oi a. differential field extension E of K(x, y) is called hyper- 
exponential (over K(j;, y)) if 

e K X, y) and € K{x, y). 

When /i G E is a hyperexponential term and ri , r2 G K(a;, y) are such that (D^h) /h = 
ri and (Dy h )/h = r2, then D^Dyh = DyD^h implies Dyri = Dxr2- Conversely, 
Christopheij ( 1999( ) has shown for algebraically closed ground fields K that for any 



two rational functions ri,r2 G K(a;, y) with Dyri — Dxr2 there exist a/h G K(a;,?/), 
Co, . . . , CL G K[x, y] and ei, . . . , G K with 

DxCq , n ("■\ , J DyCa (a\ ^ DyCi 

ri = +Dx[-r] + y and r2 = + Dy[- ] + } ei—^. 

Co \h) ci Co "Khl j-^ ci 



Together with Theorem 2 of iBronstein et al.l (|2005f ). it follows that there exists a dif- 



ferential field extension E of K(x,y) and an element ft, G E with {Dxh)/h = ri and 
{Dyh)/h — r2 which we can write in the form 



/i = coexp(|j J|c|^ 



where a G K[a;,2/], 5, co,...,cl G K[a;,y] \ {0}, ei,...,ei G K, and the expressions 
exp(a/6) and c"^' refer to elements of E on which Dx and Dy act as suggested by the 
notation. We assume from now on that hyperexponential terms are always given in this 
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form, and we use the letters a, 6, cq, . . . , cl, ei, . . . , consistently throughout with the 
meaning they have here. 



Example 2. h ^ cxp{x^y)y/x~~2y is a hyperexponential term. We have 
D^h 1 + Ax'^y - 8xy^ „ 1 , , 

^^'^^^^^^^~^eKix,y). 
n X — Zy X — Zy 

For this term, we can take cq = 1, a x^y, 6 = 1, ci = a; — 2y, ei — ^. 

We may adopt the additional condition (without loss of generality) that the q {£ > 0) 
are square free and pairwise coprimc, and that ^ N for all ^ > 0. The estimates derived 
below do not depend on these additional conditions, but will typically not be sharp when 
they arc not fulfilled. For simplicity, we will exclude throughout some trivial special 
cases by assuming that all are nonzero and that maxjdegj. a, deg^. b} + deg^, 
and maxjdegy a, deg^ b} + J2e=i '^^Sy nonzero. These latter two conditions encode 

the requirement that h is neither independent of x nor independent of y, nor simply a 
polynomial. 

Applied to the hyperexponential term h, the method of creative telescoping consists 
of finding, by whatever means, polynomials ■ • ■ iPr S IK[a:], not all zero, and a hyper- 
exponential term Q such that 

Poh+ piDxh + ■■■ +prD'^h = DyQ. 

An equation of this form is called a creative telescoping relation for /i, the differential 
operator P po + piD^ + ■ ■ ■ + prD^. appearing on the left is called the telescoper and 
Q is called the certificate of the relation. The telescoper is required to be nonzero and 
free of y, but the certificate may be zero or it may involve both x and y. When pr 7^ 0, 
the number r is called the order of P, and d := max^^Q deg^pi is called its degree. 

To motivate the form of a creative telescoping relation, assume that h ~ h{x, y) can be 
interpreted as an actual function in x and y and consider the integral f{x) = jj^ h{x, y)dy. 
Then integrating both sides of a creative telescoping relation implies that / satisfies the 
inhomogeneous differential equation 

Po{x)f{x) +pi{x)D,f{x) + ■■■ +pr{x)Dlf{x) = [g(x,2/)]JJ^^. 

In the frequent situation that the inhomogeneous part happens to evaluate to zero, this 
means that the telescoper of h annihilates the integral /. 

Example 3. A creative telescoping relation for h ~ exp(a;^?/)-y/a; — 2y is 

(3x^ - 6)/i - 2xD^h ^ Dy{{3x - 4?/)/i). 

It consists of the telescoper P = {3x^ — 6) — 2xD^ and the certificate Q = (3a; — 4j/)/i. For 
the definite integral f{x) :— exp(x'^y)y/x — 2ydy, we obtain the differential equation 

{3x^ - 6) fix) ~ 2xD,f{x) = 0. 
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Creative telescoping relations for hypere xponential terms can be found with the al- 
gorithm of Almkvist and Zeilberger ( 199Cll ). which relies on a continuous analogue of 



Gosper's summation algorithm. A more direct approach was considered by Apagodu 
(alias Mohammed) and Zeilberger ( 2005 : 2006[ ). After making a suitable choice for the 



denominator of Q, they fix an order r and a degree s for the numerator of Qj rnake an 
ansatz with undetermined coefficients, and obtain a linear system by comparing coef- 
ficients. Appropriate choices of r and s ensure that this linear system has a nontrivial 
solution, and also lead to a sharp bound on the order r of the telescoper. 

Let us illustrate this reasoning for the case where the integrand is a rational function 
h = u/v £ K(a;, y) with deg^, u < deg^ v and v irreducible. Fix some r. Then we have to 
find poi ■ ■ ■ ,Pr S K(a;) and a rational function Q € K(x, y) with 

Poh + piDxh + ■■■ + PrD^h = DyQ. 

A reasonable choice for Q is Q = (X]i=o liV'') ^ where s = deg^^ u + (?■ — 1) deg^^ v and 
qo, . . . ,qs are unknowns, because with this choice, both sides of the equation arc equal 
to a rational function with the same denominator v'^'^^ and numerators of degree at 
most degj, u + r degy u in y in which the unknowns pi and qj appear linearly. Comparing 
coefficients with respect to y on both sides leads to a homogeneous linear system of at 
most 1 -|- degy M -I- rdegyV equations with (r -I- 1) -I- (s + 1) unknowns and coefficients 
in K(a;). This system will have a nontrivial solution if r is chosen such that 



(?' + 1) + (s + 1) > dcgy u + r dcgy v + 1 r > deg 



All these solutions must lead to a nonzero telescoper P because any nontrivial solution 
with P ~ would have a nonzero certificate Q with DyQ = 0, and this is impossible 
because s was chosen such that the numerator of Q has a strictly lower degree than its 
denominator. 

We have thus shown the existence of telescopers of any order r > degy v. This is a 
good bound, but it does not provide any estimate on their degrees d. We will next derive 
inequalities involving both r and d by constructing linear systems with coefficients in K 
rather than in K(a;). 



3. Shaping the Ansatz 

Let /i be a hyperexponential term and consider an ansatz of the form 

r di / -^1 ^2 

p = E Ep^^^^'^- Q = (E E ^.^^'y' 

i;=0 j=0 ^ i=o j=o 

for a telescoper P and a certificate Q. The plan is to find a good choice for the parameters 
r, si, S2, w, do, . . . , dr- The only restriction we have is that the linear system obtained from 
equating all the coefficients in the numerator of the rational function {Ph — DyQ)/h to 
zero should have a solution in which not all the p^.j are zero. The remaining freedom can 
be used to shape the ansatz such as to keep d := max^^Q di small. 

As a sufficient condition for the existence of a solution, we will require that the number 
of terms x^y^ in the numerator of the rational function [Ph — DyQ)/h (i.e., the number 
of equations) should be less than J2l=ai^i + 1) + (*i + 1)(*2 + 1) (i-e., the number of 
variables pij and qi,j). As shown in the following example, this condition is really just 
sufficient, but not necessary. 
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Example 4. Let h = u/v be the rational function from the introduction. With r = 3, 

do = di = d2 = = d = 54, and Q = {J2i=o'l2j=o1iJ^'"y^) /""^^ comparing the 
coefficients of the numerator of {Ph — DyQ)/h to zero gives a linear system with 787 
variables and 792 equations. This system has a nonzero solution although 792 > 787. 

This phenomenon is not only an unlikely coincidence in this particular example, but it 
happens systematically when the parameters of the ansatz are not well chosen. Estimates 
which are only based on balancing the number of variables and the number of equations 
will then overshoot. It is therefore preferable to shape the ansatz for P and Q in such a 
way that the linear system originating from it will have a nullspace whose dimension is 
exactly the difference between the number of equations and the number of variables (or 
if there are more equations than variables). 

The goal of this section is to describe our choice for the ansatz of telescoper and 
certificate. The form of the ansatz for the telescoper is given in Section 3.1, the certificate 
is discussed in Section 3.2. In the beginning, we collect some facts about the rational 
functions which are used later for calculating how many equations a particular 

ansatz induces. The following notational conventions will be used throughout. 

Notation 5. • \CzP and deg^p refer to the leading coefficient and the degree of the 
polynomial p with respect to the variable z, respectively. For the zero polynomial, we 
define deg^ :~ — oo and Ic^ := 0. 

• p* refers to the square free part of the polynomial p with respect to all its variables, 
e.g., {{x+l)^{y+3fy = (a;+l)(y+3). Note that p* is only unique up to multiplication 
by elements from K\{0}, but that for any choice oip*, the degrees deg^.^* and deg^p* 
are uniquely determined and we have that p*{Dxp)/p is a polynomial in x and y. These 
are the only properties we will use. 

• := z{z- 1)(2- 2) ■■■{z-n + 1) and 2" z{z + l){z + 2) ■■■{z + n~l) denote the 
falling and rising factorials, respectively. For n < we define z— z" := 1. 

• If 2 is a real number, then 2^ := max{0, 2}. 

• If 2 is a real number, then [2J := max{a; G Z : a; < 2}, [2] := minja; G Z : .t > 2}, 
and [2] := [2 + denotes the nearest integer to 2. 

• If <I> is a formula then |$] denotes the Iverson bracket, which evaluates to 1 if $ is 
true and to if $ is false, e.g., 2+ = [2 > 0]2; 5i^j = |i = j], etc. 

Lemma 6. Let ft, be a hyperexponential term and z > 0. 
(1) If deg^. a > dcgx b, then 

DIM _ Nj 

h ~co(66*ntiC.)' 
for some polynomial Ni € K[a;,2/] with 

L 

deg^ = deg^ co + i (deg^ a + deg^ b* +^ deg^ Q - , 

L 

degy N^ < degy co + i (max{degj^ a, degy b} + degy 6* + ^ deg^ , 

L . 

Icx Ni = {Icx Co) (\cx ab* (deg^ a - deg^ bf. 

e=i 
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(2) If degj. a < deg^. 6, then 

Dih N, 



h co{bb*ULiCey 
for some polynomial Ni g K[x,y] with 

L 

deg^ N, = dcg^ Co + i (dcg^ b + deg^ + X! '^^'^x Cf - - |w G N A i > 

L 

degy N, < dcgy Co + i (ma^{dcgy a, degy b} + degy &* + ^ deg^ 

L 



ce , 



Ic.. N, 



{\Cx co)(lc^ 66* n Cf) w- if w ^ N or i < w; 



(Ic. A^c.+i)(lcx66* n Cf)'"'"+'^(-J- if w e M and ^ > 

L 

where w := deg^, co + e« deg^ ce and, if a; S N, 

L 

S := deg^ Co + (w + 1) (^dcg^ 6 + dcg^ ^* + X! '^^^^ ce - ij - deg^ A'^^+i > 1. 

e=i 

Proof. All claims are proved by induction on i. For i = 0, there is nothing to show in 
any of the cases. The calculations for the induction step i ^ i + 1 are as follows. 
(1) Let V 66* Jlfci and write for the claimed value of deg^ Ni. Then 



^coexp(^)nc?j/(coexp(^)n 



_ {D^N,)v-iN,D^v ^ {D,a)b* ~ab*{D,b)/b ^ N, ^ D^ct 
cow*+-'^ cof* 66* cow* ^-^ Ce 

{D.,N,)v ~ iN.D^v + N,{Y{ ci) {{D^a)b* - ab*^) + N,v ee^ 
e=i e=i 

cow'+-'^ 

Since deg^, a > dcg^ b by assumption, we have 

L 



deg, ({D,N,)v - iN,D^v + N,v ^ 



DxCe 

ee 



Ce 

< deg^ + deg^ v - 1 = nii + deg^ 6* + deg^. ^ + ^ deg^ C£ - 1 

L 1=1 

< mi+ deg^ a + deg^ ^* + X] '^'^^^ Q - 1 to^+i . 



e=i 

Furthermore, because of 



(£',.0)6* ~ab*^ = {\Cx a)(lc, ^*)(deg,, a - deg^ 6)a;'^°s- b'-i 



we have 



1=1 1=1 

This completes the proof that {Dl^^h)/h has the denominator as claimed and that 
its numerator has degree and leading coefficient with respect to x as claimed. The 
remaining degree bound with respect to y follows from 

^ n s 



ci 

L 



< degy co+i (degy 6*+ max{degy a, degy &}+^ deg^^ + deg^ 6*+ deg^^ '^^^v 



bounds degy Ni bounds dcg^ v 

L 

< degy Co + (i + 1) (degy b* + max{degy a, deg^ b} + ^ degy a 



and 

degy(^N,(^{D,a)b* - ab*^)]^^) 

e=i 

L L 

< degy CQ+i (degy &*+ max{degj^ a, deg^ b}+^ degy cgj + dcg^ &*+ dcg^ a+ ^ dcg^ a 



1=1 



bounds degy Ni bounds dcg^ of the other faetors 

L 

1=1 

L 



1^ 

< degy Co + (i + 1) (degy b* + max{dcgy a, dcg^ &} + ^ dcg^ < 



(2) Again, let v := bb* Y[e=i '^^ ^^'^ write rrii for the claimed value of deg^ Ni. Then, 
like in part 1, 



^,+1 , {D.:,N,)v - tN,D.,v + N,{f[ ce) {{D,a)b* - ab*^) + N,v ^ 

" _ fcl i=l 

h cow*+-'^ 

First consider the case a; ^ N or i < w. 

Since deg^. a < deg^ b by assumption, and because of 

{D,a)b* -ab*^^ (Ic, a)(lc, 6*)(deg, a - deg, b)x''''s^ "*-! + ..., 



we now have 



deg,[N,(^iD,a)b* ~ ab*^)l[, 

e=i 

< mi + deg,j. b + deg^ ^* ^ 1 + X! '^^^'■^ ^ mi+i- 



=1 

L 
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Note that this estimate is also strict when deg^ a = deg^. b because the coefficient 
of a;dog,a+dog,b*-i {D^a)b* - ab*{D^b)/b contains the factor deg^ a - deg^ 6, 
which vanishes in this case. 

Next, using the induction hypothesis, we have 

{D,Ni)v - iN,D,v + N,v ee- "'"^ 



1=1 ^ 



L 

Wi+dog^ii-l 



{\Cx Ni){\Cx v) [Acg^ Ni - i dcg^ i; + ^ deg^ q j x 

1=1 

L . L 

0-Cx Co) (ic. bb* n w- (Icj. v) (nii - i deg^ v + ^e^ dcg. 

i=i e=i 
{ICa; Co) (\Cx bb* Y[ ce) '^^(dcg^ cq + ^ deg^ a - i) a;™'+i + 



E=l 
L 



(lc,co)(lc,66*[]c,)'^'c^i±i; 



e=i 

Since w^ii 7^ when uj^'N 01 i + l<uj, this completes the proof that {Dl^^h)/h 
has the denominator as claimed and that its numerator has degree and leading 
coefficient with respect to x as claimed. The degree bounds with respect to y are 
shown exactly as in part 1. 

Now consider the case where w S N and i > uj. In this case, we start the induction 
at I = cj + 1. The induction base follows from the calculations carried out above for 
i > oj, the fact w^i^ii = 0, and the definition of S. (Note that w^itl = also implies 
(5 > 1.) For the induction step i i— > z + 1, we have, similar as before, 

deg^(^Ni({Dxa)b* - a^*-^) II c^) < "^i+i 

and 

1=1 

L 

iTii+dcg^ v — 1 



,Ct\X 



(ICa; Ni){\Cx v) {rni - i deg^ w + ^ deg^ i 

e=i 

L 

(Icx N,)(lca: v) (deg^ Co + i(deg^ v - I) - 6 - i dcg^ v + ^ee deg^ ci^l 

1=1 

(Ic. iV^+i) (ic. bb* Y[ q) '"^"^'^ i-S - l)!d^(lc, v){uj-S- i)x"'+i + • • • 



e=i 

L 



= (Ic. iV.+i) (ic, bb* n cef^' (-<5 - + . . . 

e=i 

Because of 5 > 0, the factor {—S — I'^lzl—Hl is nonzero for all i> uj. □ 
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Example 7. The case when ft, = w/w is a rational function is covered by part 2 of the 
lemma. For example, for h = {2x^ — 3x* + 5)/ {3x^ — Ax + 8) we can take cq = 2x^ — 3x^ + 5, 
a ~ 0, b — 1, L = 1, ci = 3x'^ — 4a; + 8, ei = —1. Direct calculation of the derivatives 
gives 



i 





1 


2 


3 


4 


5 


6 


Aeg,coe,{Dlh)/h 


5 


7 


9 


8 


10 


12 


14 


lc^coc\iDlh)/h 


2 


12 


36 


1512 


-18144 


272160 


-4898880 



The lemma makes no statement about the degree or leading coefficient in the case i = 
cj + 1 = 3, but knowing these, it correctly predicts all the other data in the table. In this 
example, we have 6 = 3 = uj + 1. This is not a coincidence, as we shall show next. 

Lemma 8. Let h he a, hyperexponential term with deg^ a < deg^ b, and let tu and S be 
as in Lemma 6.(2), S N. Then S > u! + 1. 

Proof. Rewrite h = cq exp{^) Yii^i c^' = cq exp{ji) Y[g^i c^' with cq = x^ , ce = cg 
{e = 1,...,L), ee = ei {i = 1,...,L), cl+i = cq, cl+i = 1, cl+2 = x, eL+2 = -i^- 
The rational functions {Dl.h)/h arc of course independent of the representation of ft, 
but the representations of these rational functions which are given in Lemma 6 are not. 
The representation obtained for the new representation of ft is obtained from the original 
representation by multiplying numerator and denominator by x^~^^c]^^ . Observe that this 
modification does not influence the values for oj or 5. It is therefore sufficient to prove the 
claim for terms of the form h = x'^h, where ft is some hyperexponential term for which 
the value of to is zero. We do so by induction on oj. For co = 0, we have S > 1 = uj + 1 
already by Lemma 6.(2). Now assume that w > is such that for x'^h the degree drop 
5 is w + 1 or more. Then for ft = x"+^ft = x{x'^h) we have D^h = x'^h + xDx{x"h), 
Dlh = 2Da:{x'^h) + xDl{x"h), and so on, ah the way down to 

i:>^+2ft = {uj + 2)D'^+\x'^h) + xD%+'^{x'^h) 
^ + 2)^^x-ft + x^^x-ft 

_ (u; + 2)N^+iv + xN^+2 j, ...s 

where N^^+i and are as in Lemma 6 and v refers to the denominator stated there. 

If 5 denotes the degree drop for ft, then this calculation implies ^ > ^. By induction 
hypothesis, we have J > w + 1. If in fact J > w + 2, then we arc done. Otherwise, if 
(5 = cj + 1, then 

Icx = {-5 - 1) Icj; N^+i Icj; w = -(w + 2) \cx N^+i Ic^; v 

by Lemma 6, so the leading terms of the two polynomials in the numerator of (1) cancel, 
and therefore 5 > uj + 1 also in this case. □ 

Experiments suggest that the bound in Lemma 8 is tight in the sense that we have 
5 = u -\- 1 for almost all hyperexponential terms ft. But there do exist situations with 
5 > w + 1. For example, it can be shown that for h = cq exp(a/6) with dcg^. b — dcg^. a > 
deg^ Co = w we have S > deg^, b — deg^ a. 
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Also Lemma 6 is not necessarily sharp for degenerate choices of h. In particular, we do 
not claim that the numerators and denominators stated in Lemma 6 are coprime. It may 
be possible to carry out a finer analysis by considering the square free decomposition of cq, 
or by taking into account possible common factors between b and the ci, or by handling 
the Ci which do not involve x separately. For our purpose, we believe that the statements 
given above form a reasonable compromise between sharpness of the statements and 
readability of the derivation. 

Several aspects of the formulas in Lemma 6 are important. One of them is that the 
denominators corresponding to lower derivatives divide those corresponding to higher 
derivatives. This has the consequence that when the linear combination Ph is brought 
on a common denominator, the degree of the numerator will not grow drastically. In a 
sense, this fact is the main reason why creative telescoping works at all. Our next step 
is to bring the formulas from Lemma 6 on a common denominator. 

Lemma 9. Let /i be a hyperexponential term and r,i E Wi with r > i > 0. 

(1) If dcgj. a > degj. 6, then 

^ co{bb*llLi^^Y 

for some N^^i G ]K[x,y] with 

L 

deg^ N„ = deg^ cq + r (dcg^ b* + deg^ ^ + X! ^'^Sx Cf) + i (dcg^ a ~ deg^ b-1 

1=1 

L 

dcgy Nra < dcgy Co + r(^dcgy b* + max{dcgy a, dcg^ b} + '^'^^v '^^) ' 

£=1 

L ^ 

Ic^ Nr^i = (Icx co)(lCa; a)'(lcx by~'(\Ca: b*y(\Cx Yl (dcg^ a - dcg^ b)\ 

1=1 

(2) If degj. a < deg^. b, then 



h co{bb*ULiC£Y 

for some Nr.i £ ]K[x,y] with 

L 

deg^ Nr^i = degj. cq + r(^dcg^ b* + dcg^ ^ + '^'^S-^ - « - € IN A i > uj\5, 

1=1 

L 

degy Nr,i < degy cq + r(^dcgy b* + max{degy a, deg^ 6} + ^ dcg^^ c^) , 

(=1 

{\Cx co)(lc^ 66*)'' (ic^ n ceYuj- if w ^ N or i < 

\c N - i ^=1 

(lc:r A^c.+i)(lc:.66* 11 ) '"^""^^^-(5 - 1) '~ ^"+^^ if w G N and i > w, 
e=i 

where w, S, and Ic^ A'cj+i are as in Lemma 6.(2). 

Proof. Both parts follow directly from the respective parts of Lemma 6 by multiplying 
numerator and denominator of the representations stated there by (bb* Y[e=i ^f)^ ^ ■ '-' 
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Since we will be frequently referring to the quantities in this lemma, it seems convenient 
to adopt the following definition. 



Definition 10. For a hyperexponential term h, let 

L 

a = deg^ b* + deg^ b + ^ deg^ cg, P = deg^ a - deg^ 6-1, 

1=1 

L L 

1 = degy b* + max{degj^ a, deg^^ ^} + '^'^^v ^ ^ + X! c<!- 

f=l £=1 

If deg^ a < deg^ b and e N, wc further let 5 be any integer with 

/ / TT- D'^+^h\ 
Lj + l<S< deg^ Co + (cj + l)(a - 1) - deg,, ( co[bb* [[ctj — j . 

Otherwise, if deg^. a > deg^ 6 or a; ^ N, let (5 = 0. Finally, we define the following flags: 

L 

03 = I f e K(.t) a V £ : (degy q = V G Z) A deg^ co > E deg^ a] . 

1=1 

Note that none of these parameters depends on r or i. The flags 0fc (fc = 1, 2, 3) are 
in {0, 1}, w belongs to K, /3 belongs to N U {—1}, and all other parameters are positive 
integers. The best value for 8 is the right bound of the specified range, but since this 
value cannot be directly read of from the input, we do not insist that b be equal to this 
value, but we allow 8 to be any number between the bound from Lemma 8 and the true 
degree drop. The flags 0i and will be used below in the ansatz for the telescoper, 03 
will play a role afterwards in the ansatz for the certificate. 

In terms of the parameters defined in Definition 10, the degree bounds of Lemma 9 
simplify to 

deg^ Nrs < deg^ co + ar + max{/3, — l}i — [oj G N A i > uj]S, 
degy Nr,i < degy cq + -fr. 

3.1. The Ansatz for the Telescoper 

Lemma 9 suggests reasonable choices for the degrees di in the ansatz for P. In partic- 
ular, our choice is based on the following features of the formulas in Lemma 9. 
• The degree of the numerator in {Dl.h)/h varies with i. A good choice for the degrees di 
will compensate for this variation, taking higher values for di when the numerator 
of {Dlh)/h has low degree, and vice versa. This is t he key idea of the Verbaeten 



completion ( Verbaeten . 1974 , 1976 : IWegschaider , 1997). 



• The leading coefficients of N^^i {i > 0) are polynomials in y, but in case 2, most of 
them are K-multiples of each other. When a and b are such that (Ic^; a)/{\cx b) € K, 
then this is also true in case 1. We will use this fact for eliminating several equations 
at the cost of a single variable. 

Before describing the ansatz for P in full generality, we motivate the construction by 
an example. 
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Fig. 2. The ansatz for P discussed in Example 11 

Example 11. Suppose that h is hyperexponential with Ic^; a = Ic^b, /3 = 1 (case 1 of 
Lemma 9), and deg^, cq = 0. 

Let r = 5 and d = 7. We want to choose d,; such that max^^Q di = 7 and the ansatz 

5 di 

leads to "many" variables but only "few" equations. The choice with most variables is 
clearly to set di — d = 7 for all i. But this ansatz leads to quite many equations. Each 
term D\. contributes to the common numerator a polynomial N^j whose degree in x 
is 5a + i + j and whose degree in y is at most 57. Because of the term x^Z?^, we must 
expect up to (5a + 13) (57 + 1) terms in the numerator. This is the expected number of 
equations in the linear system resulting from coefficient comparison. 

If we remove the term x'^ from the ansatz, i.e., if we choose do = ■ ■ ■ = d^ = 7 , 
d^ = 6, then the number of equations drops to (5a + 12) (57 + 1) because all terms x^ D\. 
other than x'^ contribute only polynomials x^ N^^ of lower degree. We save 57 + 1 
equations at the cost of removing a single variable. Removing also the terms x'^D'^ and 
x^D^ lowers the number of equations further to (5a+ 11)(57+ 1), and in general, for any 

< u) < 5, choosing di = 7 — {w + i — 5)+ (i = 0, . . . , 5) leads to (5a + 13 — ti;)(57 + 1) 
equations. The number of variables is (5 + 1)(7 + 1) — J2k=i fc = 48 — ^w{w + 1). 

If w > 1, we can introduce it; — 1 new variables by exploiting the second feature of the 
formulas in Lemma 9 as follows. Consider the choice w = 3, i.e., the terms x^ D\ with 

1 + j > 10 have been removed from the ansatz. We reintroduce the terms x'^D^, x^D'^, 
x^Dl by adding 

P3,7((deg, a - deg, bfx'^Dl - x'^DI) +P4,6((deg, a - deg, h)x'^Dt - x'^DI) 

to the ansatz, getting back the two variables ^3^7 and p^ Q but no new equations, because, 
according to Lemma 9.(1), the assumption Ic^; a ~ Ic^ b implies 

(deg^ a - deg^ 6)^ Ic^: A^5,3 = Ic^; iVs^s and (deg^ a - deg^ b) Ic^ A^5,4 = Ic^; iVs^s- 

The final ansatz is depicted in Figure 2. A bullet at (i,j) represents a variable pi_j in 
the ansatz. White bullets correspond to the reintroduced variables p^j and pi^e which 
do not affect the number of equations. 
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The general form of our ansatz for the telescoper is given in the following lemma. The 
first case is like in the example above when /3 > 0. For /3 = 0, no degree compensation 
is possible because all Nr^i have the same degree. But if {Icx a) / {ICx b) € K, it is still 
possible to save some equations by exploiting the linear dependence among the leading 
terms. In the second case, there is always a degree compensation possible, but unlike 
in the example above, terms are removed for indices i close to zero rather than close 
to r. When a; G IN, we provide an alternative ansatz which takes the degree drop 6 
into account. Common to all cases are the two basic principles of choosing di such as 
to compensate for the different degrees of the Nr^i in Lemma 9, and of installing some 
additional variables by exploiting the knowledge about the leading terms of the Nr.i- For 
the size of the cutoff, we use a new integer parameter w, whose optimal value will be 
determined later. 

Lemma 12. Let /i be a hypcrexponcntial term, r > 1, d > 0. 

(1) Suppose that dcg^. a > deg^. 6. Let < w < min{r, d//3} {w := if /3 = 0), 
di := d — /3{w + i — r)+ — 4>2 (* = 0, . . . , r), and 

1=0 j=0 

r-1 
i—r — w-\-l 

i=0 

Let N ^ co(6*&nfci CiY{Ph)/h. Then 

deg^ N < deg^ Cq + + (a + I3)r - Pw ~ 4>2 and deg^ N < deg^ cq + 7r. 

(2) Suppose that deg^, a < dcg^ b. Let < w < min{(i+ 1, r + 1}. Let di := d— (w — i)+ 
(i = 0, . . . , r), and 

r di w — 1 

i=0 j=0 i=l 

Let N = co(6*6nLi ciY{Ph)/h. Then 

deg^ N < deg^, cq + d + ar — w and deg^ N < degy cq + jr. 

(2') Suppose that deg^. a < deg^. b and cj e IN. Let uj < w < min{d — (5 + 1, r + 1}. Let 
di := d — [w — z)"*" — \i< uj\8 (i = 0, . . . , r), and 

r di UJ 

p = E E^^.^^'^^ + E^^.'^.+i {^"'^'Di - 

w — 1 

(See Figure 3 for an illustration of the shape of P in this case.) 
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Let N = co{b*bYl^^^ceY{Ph)/h. Then 

deg^ N < deg^ CQ+d + ar — w — S and deg^, < dcg.j, cq + jr. 

Proof. (1) We apply Lemma 9.(1) to each term in the ansatz for P. The claim about 
degy N follows directly from the bound on deg^, Nr.i there. For the bound on deg^ N, 
first observe that 

deg^ x^NrA < di + deg^ cq + ar + /3i 

= deg^ Co + d + ar + fii — I3{w + i — r)+ — (j)2 

= degj. Co + d + ar + I3{i — max{w + i — 0}) — 4>2 

< deg^ Co + d + ar + I3{r — w) — (j)2 

for all i,j with < i < r and < j < di. This settles the terms coming from the 
double sum. For the terms in the first single sum, which only appears when (3^0, 
we have 

deg^ x'^^+'^Nr.^ = deg^ x'^-+^Nr,r 
= deg^. Co + d + ar + /?(r — w) + 1 — 02 

and 

(fef(/3 + l))'lc.iV,,, = lc,iV,,, 
for i~r — w + l,...,r — 1. This implies 

deg. ((ilTf (/3 + IjYx'^'+'Nr^^ - x''^+^Nr,r) < deg, co + d + ar + I3{r - w) - 02, 

as desired. The argument for the second single sum, which only appears when f3 ~ 0, 
is analogous. 

(2) Now we use Lemma 9.(2). Again, the claim about deg^ N follows immediately. For 
the bound on deg^, A^, first observe that 

deg^ x^ Nr^t < di + deg^ co + ar ~ i 

= deg^ Co + d + ar — i — {w — i)'^ 
< deg^ Co + d + ar — w. 

This settles the terms in the double sum. For the terms in the single sum. we have 

deg^^ a;'*'+^iVr,i = deg^^ a;*'+^iVr,o = deg^ co + d + ar-w + 1 
and ICa; Nr^i = Ic^: Lo-Nrfi for i = l,...,u'— 1, and therefore 

deg^ [x'^^+^Nr,^ - uj^x'^°+^Nr,Q^ < dcg^ co + d + ar-w. 

(2') In this case, the terms in the double sum contribute polynomials of degree 

degj, x^Nrs < di + deg^ cq + ar — i — [i > uj]5 

= deg^ Co + d + ar — i ~ {w — i)^ — [i < — p > uj]S 
< deg^, Co + d + ar — w — 5. 

For the terms in the first single sum, we have 

deg^ x'^^+^N,.,, = deg^ x^+^iV^^o = deg^ co+d + ar-w + l- S 
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I I I 

L0 = 2 r = 11 

Fig. 3. The ansatz for P in case 2' of Lemma 12 

and Icx Nr^i ~ Ic^. uj-Nrfl for i = 1, . . . , w, and therefore 

dcg^ (^x'^^+^Nra - w-x*+i7Vr,o) < deg^ co + d + ar - w - (5. 
Similarly, for the terms in the second single sum, we have 

deg^ a;*+^7V,.,, = deg^ a^'^'+'+^A^r.^^+i < deg^ cq + + ar - w + 1 - (5. 

If the inequality is strict, we are done. Otherwise, 8 is maximal and we have 
ICa: A^r,i = lc2;(— (5 — l) '~^"''""'"'' iVr,tj+i « = w + 2, . . . , — 1 , and therefore 

deg^l^x^^'+^TV,,, - (-(5 - i) '-("+i) ^rf^+i+ijv^^^^^^ < deg, cq + d + ar - - 5, 

and we are also done. 

□ 

Lemma 12 makes a statement on the number of equations to be expected when the 
ansatz for P is made in the form as indicated. This number of equations is equal to the 
number of terms x^y^ in N ^ and this number is bounded by (deg^ iV + l)(degy + 1), 
for which upper bounds are stated in the lemma. We also need to count the number of 
variables j- . This number is easily obtained from the sum expressions given for P in the 
various cases by replacing all the summand expressions by 1. After some straightforward 
and elementary simplifications which we do not want to reproduce here, the statistics 
are as follows. 

• In case 1, the number of variables is 

(r + \){d + 1) - \p,w{w + 1) + - 1)+ - (/.2. 

• In case 2, the number of variables is 

(r + l)(d + 1) - \w{w + 1) + (w - 1)+. 

• In case 2', the number of variables is 

(r + l)(d + 1) - \w{w + 1) - + 1) + cj + (w - cj - 2)+. 
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These are only the variables coming from the degrees of freedom in the telescoper P. 
Wc will next discuss the ansatz for the certificate Q, which will bring many additional 
variables, but, by a careful construction, no additional equations. 

3.2. The Ansatz for the Certificate 

The design of the ansatz for the certificate is much simpler. Here, the goal is to set 
up Q in such a way that {DyQ)/h has the same denominator and the same numerator 
degrees in x and y as {Ph) / h docs (in order to not create more equations than necessary) , 
and that {DyQ)/h cannot become zero (in order to enforce that P 7^ in every solution 
wc find). 

A direct calculation like in the proof of Lemma 6 confirms that the first requirement 
is satisfied by choosing 



co{bb* n^^i Cf)' 



with 



deg^ Co + d + {a + /3){r ~ 1) — jSw — (j)2 — I in case 1 of Lemma 12; 
si = degj. cq + d + a{r — 1) — w in case 2 of Lemma 12; 

eg^ Co + d + a{r — 1) — w — S in case 2' of Lemma 12 

and 

S2 = degy Co + j{r - 1) + 1 in all cases. 

This ansatz provides (si + l)(s2 + 1) variables. To ensure that DyQ ^ for every choice 
of (7ij, observe that DyQ = can only happen if /i is a rational function with respect 
to y, meaning a,b G K[x] and q G K[a;] for all £ with ei ^ %. In this case, we have 
DyQ = if and only if the are instantiated in such a way that the resulting Q is 
free of ?/, and this can only happen if the choice of g^.j is made in such a way that the 
numerator degree in y is equal to the denominator degree in y. The denominator degree 
is 

L L 

^(r degy ce = 7(r - 1) - 77, where 77 = ^ deg^ 

£=1 1=1 
which is less than S2 = degy cq + 7(r — 1) + 1 if and only if deg^ cq + + 1 > 0. If we 
remove all the terms qijx'^y^ with j = 7(r — 1) — 77 from the ansatz, no instantiation of the 
remaining qij can turn Q into a term independent of y, so we can be sure that DyQ ^ 
in this modified setup. The number of variables in this modified ansatz is (si + l)s2- The 
flag (f)3 defined in Definition 10 is set up in such a way that we can in all cases assume 
an ansatz for Q with (si + l)(s2 + 1 — (^3) variables. The following lemma summarizes 
the two versions of the ansatz for Q. 

Lemma 13. Let ft, be a hypercxponential term. 
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(1) If max{degy a,degyb} > or deg^ q > for some £ with eg ^ Z, then for every 
si, S2 G N and every choice of qij € K where not ah qij are equal to zero we have 



i=0 ]=Q 



co{bb* rifci^) 



-h 



^0. 



(2) If degj^ a = degy 6 = and deg,^ cg = for all £ with eg ^ S, then for every Si, S2 € N 
and every choice of qij G K where not all qi j are equal to zero we have 



D„ 



[ Sl / ('■-l)7-')-l S2 

E E m^^^'V^ + E 

i=0\ i=0 j = (r-l)7-,,+ l 



Co (66* rifciCf) 



V 



where 7? = Efci deg C£. 



4. Solving the Inequalities 

As the result of the previous section, we obtain counts for the number of variables and 
the number of equations for a particular family of ansatzes which are parameterized by 
the desired order r and degree d of the telescoper, various Greek parameters introduced 
in Definition 10, which measure the input, and one additional parameter w by which 
the shape of the ansatz can be modulated. A sufficient condition for the existence of a 
solution of order (at most) r and degree (at most) d is 

#vars(r, d, w) — #eqns(r, d, w) > 0. 

For any particular choice of w from the ranges specified for the various cases in Lemma 12, 
we obtain a valid sufficient condition connecting r and d via the Greek parameters. 
Any of these conditions defines a region in which is inside the gray region from the 
introduction. To make this region as large as possible (and hence, as equal as possible to 
the gray region) , we will choose w in such a way that the left hand side, considered as a 
function in w, is maximal. 

It comes in handy that #;vars(r, d, w) — #eqns(r, d, w) is a (piecewise) quadratic poly- 
nomial with respect to w, so the optimal choice of w is easily found by equating its 
derivative with respect to w to zero and rounding the solution to the nearest integer. If 
this point is outside the range to which w is constrained, then the maximum is assumed 
at one of the two boundary points of the range. 

The following theorem, which is the main result of this article, contains the bounds 
which we obtained by applying this reasoning to the explicit expressions derived for 
#vars(r, d, w) and #eqns(r, d, w) in the previous section for the various cases to be con- 
sidered. 

Theorem 14. Let 

L 

/i = coexp(^)j|c|'' 

£=1 



19 



be a hyperexponential term and let a, /3, 7, S, uj, 01, (/)2, 03 be as in Definition 10 and set 
-0 = 7 + 03 — 2. Then a creative telescoping relation for h of order r and degree d exists 
whenever 

t9 r + 09 

r >■(/; + 1 and d> 



where •§ and (/3 arc defined as follows. 
(1) If dcg^. a > dcg^ b, let 



?9 = (a + f3){2-/ - 1 + 03) + 7 - 1, 

LP = deg^ Co + (a + /3 + 1) deg^^ cq + (7 - 2 + 03)(deg^ cq - a - ^ - 02) 
- (1 - 02)(7 - 2 + 03)+(0i + i/3(7 - 1 + 03)). 
(2) If degj. a < deg^. 6, let 

i9 = a(27 - 1 + 03) - 1, 

(f = dcg^ Co + a dcgj^ cq + (7 - 2 + 03)(dcg2, co + 1 - a) 
-i(7-2 + 03)+(7 + l + 03). 
If furthermore a; S N and 7 — 1 + 03 > cj and S = lo + 1, then if can be replaced by 
ip' ^ ip - (5(7 — 2 + 03 - w) + 1. 

Proof. (1) Suppose deg^ a > deg^ 6. According to the calculations done in the previous 
section, in this case there exists an ansatz with 

(r + l){d + 1) - ^pw{w + 1) + 0i(w - 1)+ - 02 

variables coming from the telescoper P, 

(deg^ co + d+(a + /3)(r - 1) - pw - 02) (deg^^ cq + 7('^ - 1) + 2 - 03) 

variables coming from the certificate Q, and 

(deg^ Co + d+ {a + I3)r - (3w - (f)2 + l) (degy co + 7r + l) 

equations. Therefore, a creative telescoping relation exists provided that 

(r + l){d + 1) - ^(3w{w + 1) + 0i(w - 1)+ - 02 

+ (deg^ Co + d + (a + /3)(r - 1) ~ ~ 02)(degy co + 7(?- - 1) + 2 - 03) 
- (deg^, Co + d + (a + /3)r - /^w - 02 + l)(dcgy co + 7r + 1) > 0. 

For ?' > 7 — 1 + 03 , this inequality is equivalent to 

d > (((a + /?) (27 - 1 + 03) + 7 - l)r + deg,, cq + {a + (3 + 1) degy cq 

+ (7-2 + 03)(deg,co-a-/?~02) (2) 
+ ^(3w{w - 27 + 3 - 203) - 0i(«; - 1)+) / (r - 7 + 2 - 03) . 

The choice w = proves the claim when 02 = 1 or 7 < 1 — 03. Now suppose 
that 02 = and 7 > 1 — 03. The claimed estimate is obtained for the choice 
'u; = 7 — 1 + 03 >0. We have to show that this choice is admissible, i.e., that 
1<7— l + 03< min{r, d/P}. Because of 7 > 1 — 03, the lower bound is clear, and 
r > 7 — 1 + 03 holds by assumption. To see that 7 — 1 + 03 < d//?, observe that 
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the right hand side of (2) converges to (a + /3)(27 — 1 + ^3) + 7 — 1 for r — > 00. 
Since its numerator is nonnegative (as checked by a straightforward calculation), 
it follows that this inequality implies 

d > (a + /3)(27 - 1 + 03) + 7 - 1 > ^(7 - 1 + </)3), 

as desired. 

(2) Now assume deg^, a < dcg^ b. From the counts of variables and equations in the 
ansatz described in Lemma 12.(2), we find that a creative telescoping equation 
exists provided that 

(r + l){d + 1) - iw(w + l) + {w- l) + 

+ (deg^ co + d + air -1) -w + l)(dcgj^ co + 7(7' - 1) + 2 - (^3) 
- (deg^ Co + d + ar -w + l)(deg,y cq + jr + I) > 0. 

For r > 7 — 1 + 03, this inequality is equivalent to 

d > (^(a(27 - 1 + 03) - l)r + deg^, ca + a degy cq + (7 - 2 + 03)(deg^ cq + 1 - a) 

+ (i - 7 - cl)3)w + -{w- 1)+) / (r - 7 + 2 - 03) . 

Regardless of the choice of w, the right hand side is at least Q!(27 — 1 + 03) — 1. 
Similar as before, the claimed bound follows on one hand from the choice w = 
and on the other hand, if 7 > 1 — 03, from the choice w = 7 — 1 + 03, which also in 
this case is in the required range because 1<7— 1 + 03 < a (27 — 1 + 03) — l<d 
and 7 — 1 + 03 < r. 

The second estimate is obtained from the alternative ansatz from Lemma 12.(2'). 
The inequality in this case is 

(r + l){d + 1) - \w{w + 1) - (5(cj + 1) + a; + (w - cj - 2)+ 



+ (deg^ Co + d + a{r ~ 1) 



5 + l)(deg^ CO + 7(?- - 1) + 2 - 03) 



- (deg^ CQ + d + ar-w-5 + l)(degj^ co + 7r + 1) > 0, 
which for ?' > 7 — 1 + 03 and ui = 7 — 1 + 03 is equivalent to 

(a(27-l + 03)-l)r + (^' 
r - 7 + 2 - 03 

It remains to show that the choice w = 7 — 1 + 03 is compatible with the range 
restrictions for w applicable in the present case. While the requirements < 7 — 
1 + 03 ^ + 1 are satisfied by assumption, the requirement 7— l + 03<(i — (5 + 1 
is less obvious. A sufficient condition is 

(a(27-l + 03)-l)r- 



7- 



> 7 - 2 + 03 + (5. 



It can be shown ea s ily with Collins's cylindrical algebraic decomposition algo- 
rithm ( Collind . I975I: Cavi ncss an d Johnson . 1998[) (e.g., with its implementation 
in Mathematica ( Strzebohs ki. 2000, 2006[) ) that this latter inequality follows from 
deg^co > 0, degy co > 0, a > 1, r > + > uj+1 >l,d = cj + 1, 03(03-1) = 0, 
and 



(p' = deg^ Co + a deg^^ cq + (5a; + 1 + (7 - 2 + 03)(deg^ co - a - 5(7 - 1 



S). 
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This completes the proof. 
□ 

As we do not claim that our bounds are sharp, no justification for the various choices of 
w are required in the proof. But of course, the choices were made following the reasoning 
outlined before the theorem. For example, in case 1 the main inequality is 

(r + + 1 + 02) - \Pw{w + 1) + 0i(w - 1)+ - 02 

+ (deg^ Co + d + (a + p){r - 1) 02 + l)(deg^ cq + -i{r - I) + 2 ~ 03) 

- (deg^ Co + d + (a + /3)r - /3w - 02 + l)(degy co + 7?- + 1) > 0. 
Differentiating the left hand side with respect to w gives 

-/3u;- 1/3 + ^7 + 01 +;303, 

which vanishes for ^ = 7— | + 03 + 0i//3. The unique nearest integer point is [7 — 
f + 03 + 0i//3] =7—1 + 03 when (jji/ (3 ^ 1. When 0i//3 = 1, there are two nearest 
integer points 7—1 + 03 and 7 + 03 , and since the maximum is exactly between them and 
quadratic parabolas are symmetric about their extremal points, the values at 7 — 1 + 03 
and 7 + 03 agree. In conclusion, the choice w = 7 — 1 + 03 is optimal in both cases. 

The calculations for the other cases are similar. But note that having chosen w op- 
timally does not imply that the bounds given in the Theorem 14 are tight, because the 
whole argument relies on counting variables and equations for the particular ansatz fam- 
ily introduced in Section 3, and we cannot claim that this shape is best possible. Recall 
that we aim at an ansatz for which the number of solutions of the resulting linear system 
is equal to (or at least not much larger than) the difference between number of variables 
and number of equations. One way of measuring the quality of our ansatz, and hence the 
tightness of our bounds, is to compare the region of all points (r, d) where an ansatz for 
order r and degree d actually has a solution (the "gray region" from the introduction) 
with the region of all points (r, d) for which Theorem 14 guarantees the existence of a 
solution. The following collection of examples shows that there are cases where Theo- 
rem 14 is extremely accurate as well as cases where there is a clear gap between the 
predicted shape and the actual shape of the gray region. As a reference ansatz for exper- 
imentally determining in the examples whether a specific point (r, d) belongs to the gray 
region, we checked whether the naive ansatz where d^ = di = ■ ■ ■ = dr (i.e., w = 0) as a 
solution, because every solution of some refined ansatz with w > is also a solution of 
the ansatz with w = 0. It is not guaranteed however that this ansatz covers all creative 
telescoping relations. Additional relations at points (r, d) outside of what we indicate as 
the gray region may exist. For example, when our ansatz leads to a solution (P, Q) in 
which all the polynomial coefficients of P share a nontrivial common factor / S K[a;], 
then (i-*//, Q/ f) is another relation with a telescoper of lower degree. This phenomenon 
can often be observed for the minimal order telescoper, but as we do not know of any 
efficient way of detecting it also for the nonminimal ones, we can unfortunately not take 
it into account in the figures. 

Example 15. (1) Consider the term h = wexp(w) where 

u = 7x^y^ + 8x^y^ + 9x^y + + lOx'^y^ + 2x'^y^ + 3x^y + 9x'^ 
+ 7xy^ + 4:xy'^ + 5xy + 3a; + 9y^ + + 6y + l, 
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V = 6x^y^ + 4x^y^ + x^y + 9x^ + Sx'^y^ + Sx'^y'^ + 2x^y + Sx^ 

+ 3a;y^ + Txy'^ + '^xy + %x + 5y^ + 2y^ + 7y + 6. 

We are in case 1 of Theorem 14 and have a = 0, /3 = 2, 7 = 3, = ^2 = (As = 0, 
degj. Co = degj^ cq = 3. According to the theorem, we expect creative telescoping 
relations for all {r,d) with r > 2 and d > (12r + ll)/(?' — 1). Figure 4. (a) depicts 
the curve (12r + ll)/(r — 1) together with the gray region. In this example, the 
gray region consists exactly of the integer points above the curve: the bound is as 
tight as can be. 

(2) Now consider the term h = cxp{u)/v where 

u = 4a;^y^ + 7x'^y + 9x'^ + bxy'^ + 2xy + 3x + + y + 6, 
V = 6a;^y2 + lOx'^y + + 9xy'^ + 5xy + 8x + 8y'^ + lOy + 8. 

We arc again in case 1 of the theorem and we have a = 2, ^ = 1, 7 = 4, 01 = 
(/)2 = 03 = 0, dcg^ Co = degj^ cq = 2. The estimate from Theorem 14 is now 
d > (24r— 9)/(r— 2), which is depicted together with the gray region in Figure 4.(b). 
In this case, the bound is not sharp. 

(3) Now let h be the rational function from the introduction. Then we are in case 2 
of the theorem and we have a = 3, /3 = —1, 7 = 3, w = —1, 6 = 0, (f>i = 1, 
(/)2 = Oi <^3 = Ij deg^ Co — degj, Cq ~ 2. The bound from the theorem is now 
d > (17r + 3) /{r — 2), which is shown together with the gray region in Figure 4.(c). 
The curve correctly predicts all the degrees except for the minimal order recurrence, 
where the true degree is one less than predicted. 

(4) Next, let h — u/v with 

u = Ax^y'^ + Tx^y + 9x^ + 5xy^ + 2xy + 3x + + y + 6, 

V = (6x^2/^ + lOx'^y + 6x'^ + 9xy^ + bxy + 8x + 8y^ + lOy + 8) 

X (8x^2/^ + 7x^y + 4.1;^ + 5xy^ + 3xy + 7x + 9y^ + 7y + 7) . 

This term is also covered by case 2 of the theorem, and we have a — 4, /S = —1, 
7 = 4, w = -2, (5 = -1, 01 = 1, 02 = 0, 03 = 0, deg^^co = deg^^ cq = 2. The 
estimate d > {27r + 3)/ {r — 2) from the theorem is correct but not tight, as shown 
in Figure 4.(d). 

(5) Finally, let h ~ '/u with 

u = 4x^2/6 + ^x^y^ + 2x2/ + 7x^y^ + 7x^y^ + 2x^y + 7x^ + lOxy*^ + 7xy'^ + 9xy'' 
+ Axy^ + 5xy2 + bxy + 7x + Ay^ + 3y^ + 2y^ + 8y-^ + 3y^ + 7y + 2. 

Now the alternative bound of case 2 with ip' in place of ip is applicable because we 
have w = 1 G N. The bound using (p is d > (21r — 18)/(r — 4). The first correctly 
predicted degree occurs at r = 14. In contrast, the bound d > (21r — 23)/(r — 4) 
using (/?' is tight for all r > 5 and only off by one for the minimal order r = 5. The 
situation is shown in Figure 5. On the right, we show a comparison of the sharp 
bound based on (p' (solid), the bound based on (p (dashed) and the bound which 
would be obtained by choosing w = instead of w = 7 — 1 + 03 in the proof of 
Theorem 14 (dotted). 

There are several ways of refining the ansatz for P and Q even further in order to 
achieve better estimates where ours are not sharp. Here are some ideas. 
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Fig. 4. Sizes (r, d) of creative telescoping relations together with the curve predicted by Theo- 
rem 14, for the hyperexponential terms discussed in Example 15. 

• The possibility of introducing extra variables without increasing the number of equa- 
tions (depicted by the white bullets in Figures 1 and 2) rests on the observation made 
in Lemma 9 that the leading coefficients Ic^N^.i arc K-multiples of each other, i.e., 
that these leading coefficients generate a linear subspace of K[?;] of dimension one. Ex- 
periments suggest that this observation can be generalized to the coefficients of lower 
degree as follows: If Vj C K[y] denotes the vector space generated by the coefficients 
of I'^^'Sx Af^.-i in Nr,i {i = 0,...,r), then Vb C Vi C • • • C and dim Vj < j + 1 at 
least for small j. If this is true, it would allow adding more extra variables without 
increasing the number of equations. 

• In general, comparing coefficients of the monomials x^y^ of a polynomial S to zero 
results in a linear system with (deg^. S + l)(degj^ S" -f 1) equations. But if S contains 
some factor which is free of the variables pij and qij, then canceling this factor before 
comparing coefficients results in a system with fewer equations and the same number of 
variables. While in our case, it is too much to hope for a factor which would divide S as 
a whole, it seems that at least in some cases, factors can be removed from Ic^^ S G K[y] 
or Icj, S* € K[a;]. For example, when deg^ a > deg^, 6 and deg^ a > deg^, 6, it can be 
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Fig. 5. Left: Sizes (r, d) of creative telescoping relations together with the curve predicted by 
Theorem 14, for the term discussed in Example 15.(5). Right: a detail of the figure on the left in 
a larger scale, together with the curve based on tp instead of tp' (dashed) and the curve based on 
w = Q (dotted). The correct degrees are precisely the smallest integers strictly above the solid 
curve. The two variations both overshoot for all the points in this range. 

shown that Ofci ci \ \Cx S and Jlfci ^'^v I ^^v '^^ ^'^ J2e=i i'^^&y ^'^x eg + deg^ Icj; ce) 
equations can be discarded in this case. 

We have not worked out the influence of these variations in full generahty, but only 
on some examples. It turned out that they indeed lead to tighter estimates, but the 
difference is rather small, and decays to zero for large r. At the same time, they would 
lead to much more complicated formulas. We do not know the reason for the gap in 
Examples 15.(2) and 15.(4) between the curve from Theorem 14 and the boundary of 
the gray region for r — > oo. Even though it appears more important for a bound to be 
tight for small orders than for large ones, we would be very interested in seeing a refined 
bound which closes this gap. 

It is also interesting to compare the gray regions for hyperexponential terms com- 
posed from dense random polynomials with the gray regions for hyperexponential terms 
of the same shape that originate from some specific application. According to our experi- 
ments, the shape of the gray region for a randomly chosen term h = cq exp{a/b) Y[e=i 
only depends on the number L of factors in the product, the degrees of the polyno- 
mials a, 6, Co, . . . , cl, and the exponents ei,...,eL. However, input containing sparse 
polynomials or polynomials which in some other sense have a "structure" may well have 
considerably smaller degrees. 

Example 1 6. If an,k denotes the number of HC-polynomioes with n cells and k rows 
(|Wilj . ll989l . Section 4.9), then 



,,fe 



kx y 



xy{l — x)^ 



n,k=0 



(1 - a-)4 - xy{l 



x'^y)' 



A differential equation for the generating function X^^^o '^n,nx" of the number of HC- 
polynomials with n cells and n rows can be obtained by applying creative telescoping to 
the rational function obtained from the rational function above by substituting x by y, 
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Fig. 6. Gray regions for the two terms h (light gray) and g (dark gray) from Example 16. 
Although all Greek parameters have the same values for h and g (and hence, Theorem 14 gives 
the same degree estimation curve), the actual gray regions differ significantly. 

y by x/y, and dividing the result by y. Let thus 

1 y^{^-y? ^ x(i-yf 

y (1 - yY - y|(l -y-y^ +y^ + y^^) - yf - x{l - y + xy - y"^ + y^)) ' 

Here we have cq = x{l~y)^, a 0, 5 = 1, ci = y, C2 = ((1 -y)"* -y + xy-y^ +?/^), 
ei =62 = —1. The gray region for h is shown in light gray in Figure 6. For comparison, 
the same figure contains the gray region (in dark gray) for a term g which was obtained 
from h by replacing cq and C2 by dense random polynomials with deg^. cq = 1, degj^ cq = 3, 
deg^: C2 = 2, degj^ C2 = 4, so that all the Greek parameters have precisely the same values 
for g and h. 

Theorem 14 predicts relations whenever d > ^''jZ^ (black curve), which is a good 
estimate for the generic term g but a significant overestimation for the special term h. 




5. Consequences and Applications 



Ou r theorem contains as a special case Theorem cAZ of lApagodu and Zeilberger 
(|2006[ ). which says that a (non-rational) hypercxponential term always admits a tele- 
scoper of order r = 7 + 1, but makes no statement about its degree d. Similarly, we can 
also give an estimate for the possible degrees d without paying attention to their orders r. 

Corollary 17. (1) For every hyperexponential term /i, there exists a creative telescop- 
ing relation of order r = '(/; + l = 7 + l — (/)3. 
(2) For every hyperexponential term /i, there exists a creative telescoping relation of 
degree 



(a + ^)(27 - 1 + 03) +7 if deg^ a > deg^ b; 
a(27-l + 03) if deg^ a < dcg^ 6. 

Proof. Both claims are immediate by the formulas given in Theorem 14. □ 
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In connecting order r and degree d into a single formula, Theorem 14 makes a much 
stronger statement than this corollary. Assuming for simplicity that the bounds of The- 
orem 14 are tight, we can use them to compute optimal choices for order and degree of 
the telescoper. There are various quantities which one may want to minimize. Besides 
asking for a bound on the minimal order or the minimal degree, as carried out above, 
we may ask for a choice (r, d) where the computational cost is minimal, or the total size 
S{r,d) := (r + l){d + 1) + (si + l)(deg^ cq + j{r - 1) + 2) of the output (consisting of 
telescoper and certificate), or the size T(r, d) := (r + l){d+ 1) of the output telescoper 
alone. Or, if the telescoper P is to be transformed into a recurrence for the series coef- 
ficients of its solutions, one may want to minimize the order of this recurrence, which is 
bounded by R{r,d) := r + d (see, e.g., Thm. 7.1 in lKauers and Paulell201l[) . 

For minimizing the computational cost, wc first have to fix a particular algorithm for 
computing P and Q for given h. We are not forced to follow the algorithm which is 
implicit in the analysis of Sections 3 and 4 (making an ansatz, comparing coefficients 
with respect to x and y to zero, and solving a linear system of equations over K). 
In fact, this algorithm has a rather poor performance. It is much better to do a co- 
efficient comparison with respect to y only and to solve a li near system of equations 
over ] K(x). This is also what is proposed in the original articles ( Almkvist and Zeilberger , 
19901 : Mohammed and Zcilbergeil 12005 ; Apagodu and Zeilberger , 2006 ) and what is used 



in practice (fKou ts c han , 200l l2010|). Output sensitive linear system solvers based on 



Hermite-Pade approximation (jBeckermann and Labahnl . Il994i IStoriohann and Villard 



2005: Bostan et al 



2007|) are able to determine the degree n solutions of a linear sys- 
tem over K(a;) with m variables and at most m equations using 0~(7ito^) operations 
in K. Since an ansatz over K(x) will have only r + I variables coming from the tele- 
scoper, dcgy Co + j{r — 1) — 03 -t- 2 variables coming from the certificate, and a so- 
lution of degree si with respect to x, it seems reasonable to assume that the com- 
putational cost is minimal for a choice {r,d) which minimizes the function C{r,d) := 
si(dcg^ Co + (7 + l)r - 7 - 03 + 3)3. 

Example 18. Consider a hyperexponential term h = cq exp{a/b)yfci where a, &, co, ci S 



K[a;,y] have the degrees deg^, a 
deg^ ci = 4, degj^ ci = 
a; = 4, (5 = 5, 01 = 



deg. 



a = deg^ b = degy 6=1, deg^ co = deg^^ co = 2, 
6. We are in case 2 of Theorem 14 and have a = 6, /3 = — 1, 7 = 8, 
02 = 0, 03 = 0. According to the theorem, a creative telescoping 

1 



relation exists for (r, d) with r > 7 and d > (89r - 40)/(r - 6) + 1 = (90r - 46)/(r - 6). 

On the curve d = (90r — 46)/ (r — 6), the cost function C(r, d) = (6r + d— 16) (9r — 3)'^ 
assumes its minimal value for r = 8 rather than for the minimal order r ~ 7. Finding this 
optimal value is easy: regard r temporarily as real variable and use calculus to determine 
the minimum of C(r, ^^^Zq^ ) ■ This gives a minimum point near r = 7.679. It follows that 
the minimum for r S N is either at r = 7 or at r = 8. Comparing the actual values of C 
at these two points indicates that the 8th order telescoper is about 8% cheaper than the 
7th order operator, and hence the cheapest operator of all. 

By similar calculations, we find that the output size (telescoper and certificate com- 
bined) is minimized for r = 10, the size of the telescoper alone is minimized for r = 12, 
and the order of the recurrence associated to the telescoper is minimized for r = 28. See 
Figure 7 for an illustration. 
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Fig. 7. Points (r, d) on the curve for which (a) the order, (b) the computational cost, (c) the 
size of telescoper and certificate combined, (d) the size of the telescoper only, (e) the order of 
the recurrence corresponding to the telescoper, and (f) the degree is minimal. 

For the moment, the term h considered in the above example is a bit too big to actually 
compute the creative telescoping relations of orders 7 and 8 and compare the difference of 
the timings to the predicted speedup of 8%. On smaller examples, the minimal (predicted) 
complexity is achieved for the minimal order operator. It may seem that an improvement 
by just a few percent is not really worth the effort. But in fact, the improvement gained in 
the example is just the tip of an iceberg. Asymptotically, as the input size increases, the 
speedup becomes more and mo re significant. In the next result, which is a generalization 
and a refinement of a result of Bostan et al. 1 (|2010l) . we give precise estimates. 



Corollary 19. Let h he a hyperexponential term and r = max{a, 7, deg^ Cg, degj^ Cq}. 
Let K be an increasing sublinear function with the property that degree n solutions of 
a linear system with m variables and at most m equations over K(a;) can be computed 
with nm^K{m.a,x{n,m}) operations in K. Then a creative telescoping relation of order 

r = T — I + (f>3 can be computed using 

2K(2r3)T^ + 0~(T«) 

operations in K. If r is chosen such that 

r = 1(1 + v^)t + 0(1) < 1.281T + 0(1) 

then a creative telescoping relation of order r can be computed using 

^(349 + 85VT7)k{11t^)t^ + O(t^) < 21MK{nT^)T^ + 0~(r^) 

operations in K. In particular, creative telescoping relations for hyperexponential terms 
can be computed in polynomial time. 
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Proof. First assume deg^. a > deg^ b. According to Theorem 14, there exists a creative 
telescoping relation of order r and degree d whenever r > r — 1 + (/>3 and 

a > f{r) — , 

r - T + 2 - 03 

where the term O(t^) is independent of r. A creative telescoping relation of order r and 
degree d can be computed using at most 

C(r, d) = ((r + l)r + 3 - 03)'((/S + r)r + d - /3(t + ^3) - 02 - + r){r + 1) + d) 

operations in K. The claim follows from evaluating C{r,f{r)) aX r = t — 1 + (j)y, and 
r = i(l + \/T7)r + 0(1), respectively, and replacing the arguments of k by generous 
upper bounds. 

For the case deg^ a < deg^ b, the estimates are proved analogously. Although the 
formulas for /(r) and C{r,d) are slightly different in this case, the final result turns out 
to be the same. We leave the details to the reader. □ 

The strange constant |:(1 + VT?) in Corollary 19 is chosen such as to minimize the 
multiplicative constant in the complexity bound under the simplifying assumption that 
K is constant. It was determined by first equating ^C(r, /(r)) to zero, which yielded the 
optimal choice of r as an algebraic function in r, /3, and ^3. The term |:(1 + ■\/17)t is the 
dominant term in the asymptotic expansion of this function for r — > 00. It is perhaps 
noteworthy that the choice of the constant is irrelevant for achieving a cost of 0"^(t^), 
as long as the constant is greater than I. Taking r = ut for arbitrary but fixed u > 1 
leads to the complexity bound " l^^i'^ k-t^ + 0~(r^). The choice u = j{l + \/T7) only 
minimizes the leading coefficient. Since ■|(l + vT7) ~ 1.28, the result indicates that when 
a and 7 are large and approximately equal, it appears to be most efficient to compute a 
telescoper whose order is about 30% larger than the minimum order. 

In the same way as exemplified in Corollary 19, we have also determined the choices 
for r for which some other quantities become minimal. The resu lts are given in Tabl e 1. 

As a final application, we improve some of the results given bv lBostan et al. I ()2007l) on 



differential and recurrence equations related to algebraic functions. Let m € ]K.[x,y] be 
irreducible with deg^ to > 1, and let a G K[[a;]] be such that m{x, a{x)) = 0. According to 
Proposition 2 in their paper, if P + DyQ is a creative telescoping relation for y{Dym)/m, 
then Pa = 0. Thus we can use our results about creative telescoping to derive estimates 
for differential equations for a. 

Corollary 20. Let m e K[a;,y] and a = X)i^o ^ ^[N] be as above and write 

Tx := degj. TO, Ty := degy to. Assume > and Ty > 0. Then 

(1) The series a satisfies a linear differential equation of order r ~ Ty with coefficients 
of degree 

d = 2TxTy - \Ty + T^Ty - ^Ty +Tx+3. 

(2) The series a also satisfies a linear differential equation of order r = 2tj, with coeffi- 
cients of degree 

1 r Tt + 1 ' 

d ^ AT^Ty - 2 - 3r3, - 1 - 
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(c) 
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11+5V5 4 
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(4 + 2^/5)r3 


(3 + x/5)r2 


(3 + v^)r2 


(d) 
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8r^ 


4r2 
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(e) 






2r5 


2^/2r^/2 


2t^ 


2r2 


(f) 


2r^ 




4r« 


4r^ 


2t^ 


2r2 



Table 1. Minimizing various functions on the curve of Theorem 14. The table shows the order r, 
the complexity C(r, d), the output size S{r, d) of telescoper and certificate, the output size T(r, d) 
of the telescoper only, the recurrence order R{r,d), and the degree d of the telescoper when r 
is chosen such that (a) r is minimal, (b) C(r, d) is minimal, (c) Sir, d) is minimal, (d) T{r, d) 
is minimal, (e) R{r, d) is minimal, (f) d is minimal. The parameters r and «: have the same 
meaning as in Corollary 19. The arguments of n are suppressed. Only the dominant terms of 
the asymptotic expansion for r — >■ cxa are given. In rows (e) and (f), the values for d differ only 
in the lower order terms. 

(3) The coefficient sequence (a„)^g satisfies a linear recurrence equation of order 

with polynomial coefficients of degree 

Ty-1 + ^^{8t^ - + 4)r, - 2t2 - Gr^ + 12 . 

Proof. For h — y{Dym)/m we have dcg^. cq < a ~ t^, dcgj, cq = ^ = Ty, lu < 5 < 1, 
and 03 = 1. According to Theorem 14.(2), a creative telescoping relation of order r and 
degree d exists provided that r > Ty and 

^ ^ ATxTyr + 2TxTy - Ty - ?>Ty + + 6 

- 2(r -Ty + l) ■ 

Parts 1 and 2 follow from here by setting r — Ty or r = 2tj,, respectively For part 3, 
observe first that there exists a creative telescoping relation of order r and degree d where 

r>Ty-l + iy^(8r2 - 4Ty + 4)r, - 2t^ - Gry + 12, 

d > 2T,Ty + iy^(8r2 - 4Ty + 4)t, ~ 2t^ - 6Ty + 12. 

From here the claim follows by the fact that when a power series a satisfies a linear 
differential equation of order r and degree d, then its coefficient sequence satisfies a 
linear recurrence equation of order r + d and degree r. □ 

These results are to be compared with the corresponding results of Bostan ct al. 
(degree 4r^r2 + smaller terms for part 1, order 6Ty and degree ir^Ty for part 2, and order 
and degree 2TxTy + + 1 for part 3), as well as with the conjectures about the minimal 
sizes they found experimentally {2t^ — Sr^ + 3r for part 1 when = Ty =: t and order 
and degree 2TxTy — 2 — (jx — Ty) for part 3 if Ty > 1). 
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6. Conclusion 



What is the shape of the gray region? Where docs it come from? And how can it 
be exploited? — These were the guiding questions for the work described in this article. 
As a main result, we have given in Theorem 14 a simple rational function whose graph 
passes approximately along the boundary of the gray region, in some examples more 
accurately than in others. This curve was derived from a somewhat technical analysis 
of the linear systems resulting from a specific ansatz over K. Where the curve does not 
describe the gray region accurately, these linear systems have solutions despite of having 
more equations than variables. Some possible reasons for this phenomenon were taken 
into account in the design of the ansatz, thereby improving the accuracy of the estimate 
compared to a naive approach. However, as shown in Examples 15.(2) and 15.(4), there 
seem to be further effects which sometimes cause a gap between the true degrees and our 
prediction. It would be interesting to know what these effects are, and to derive sharper 
estimates from them. Ultimately, it would be desirable to have a version of Theorem 14 
which is generically tight. 

Tight curves allow for optimizing computational cost, output sizes, and other mea- 
sures by trading order against degree. As the degree decreases when the order grows, it 
is not always optimal to compute the minimal order operator. In Example 18, we have 
illustrated how the curve of Theorem 14 can be used to calculate a priori the optimal 
orders for several interesting measures. Of course, if the curve is not tight, these predic- 
tions may not be correct, but even then, at least they provide some useful orientation. 
Tightness of the curve is also not required for deriving asymptotic bounds on the com- 
plexity. As we have shown in Corollary 19, the difference between the optimal choice 
and other choices is significant for asymptotically large input size. We believe that this 
result is not only of theoretical interest. Even if the minimal cost may be achieved for the 
minimal order in any example which is feasible with currently available hardware, it can 
be seen from Example 18 that it already starts to make a difference for inputs which are 
only slightly beyond the capability of today's computers. We therefore expect that the 
technique of trading order for degree will help to optimize the performance of efficient 
implementations of creative telescoping in the near future. 
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